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TECHNICAL MEMORANDUM 


ET LOX MODAL SURVEY ANALYSIS AND 
TEST ASSESSMENT 

I. INTRODUCTION 


The requirement that no "pogo" effects would be present during 
flight of the Space Shuttle made it mandatory that accurate dynamic 
models be available for analysis. To accomplish this task, a new 
methodology (asymmetric hydroelastic analysis) woidd be necessary. 

The Martin Marietta Corporation (MMC) developed a three-dimensional 
finite element hydroelastic capability to perform the dynamic analysis. 

The Marshall Space Flight Center (MSFC) set forth requirei.,cnts and 
performed a liquid oxygen (LO 2 ) modal survey for the LOX portion of 

the External Tank to verify the math model and to determine what, if 
any, anomalies might be present in the structure. 

The test airticle consisted of a full-scale LOg tank and flight inter- 

tank mated to a stiff support test ring fixture. The test article was 
supported in a soft spring mode by 33 airbags arranged into three 
groups (Fig. 1). 

The test conditions were selected to be representative of the flight 
environment. Four conditions were selected to represent liftoff. Solid 
Rocket Booster (SRB) separation, midrange flight, and end burn. Since 
the thrust angle of the Space .Shuttle is maintained through the c.g. , 
a canted angle of 13® was seliH.-!(>«l for the last three test conditions. 
Additional test conditions were added as the test progressed to obtain 
a better understanding of the low damping observed in the test results 
of the second bulge mode. 

The LOg modal survey begun in February 1978 and was completed 

in July 1978. The pretest analysis of the test article and test support 
structure proved to be an excellent representation of the test results of 
the tank and simulated propellant (H 2 O). Minor changes were required 

in the SRB cross beam model to match the test data. 
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II. MATH MODEL 


A. Modeling Overview 

The selection of a test article and support system was followed by 
a mathematical modeling exercise to exactly duplicate the structure and 
test supporting equipment. These models will determine the final fre- 
quencies, mode shfl|>es, and dynamic pressure. 

A new analysis technique was developed by Martin-Marietta Com- 
pany (MMC) who was responsible for the analysis and model /test data 
correlation. This analysis, an asymmetric finite element hydroelastic 
capability, was formulated by using the standard structural finite 
elements and a potential function approach with assumed incompressibifity 
for the fluid. 

A backup for this modefing capability was undertaken by tne 
Systems Dynamics Laboratory at MSEC to assure *^odeling success. A 
hydroelastic asymmetric analysis c£^abifity was developed for the MSEC 

version of NASTRAN^^\ This capability is available to NASTRAN users 
through the COSMIC Library at the University of Georgia. 


B. Structural Model 

A convenient plane of symmetry of the test article and support 
hardware allowed one- half of the test article to be modeled. The nodes 
(grid points or collocation points) are at the intersections of 13 longi- 
tudinal and 29 transverse planes. The longitudinal planes are 15® apart, 
Eigure 2. The vransverse planes were established by structural features 
such as weld lands, reinforcing rings, structural interfaces, etc. Addi- 
tional nodes were incorporated where necessary to model the SRB cross 
beam, test support ring, and suspension system. 

1. Liquid Oxygen Tank (LO 2 Tank). The LO^ tank is a three- 

part monocoque aluminum stx*ucture (Fig. 3) consisting of a half oblate 
spheriod lower dome, a cylindrical mid -section and a forward ogive 
section. A Y-shaped (Y) ring connects the cylinder /lower dome, 
supports the slosh baffles, and is the bolt ring interface for the LO^ 

tank to the intertank. A second ring. T ring, stiffens the tank, supports 
the slosh baffles, and forms the interface where the ogive is welded to 
the cylinder. A line of nodes was located at the T ring and Y ring. 

The location of the other node lines (transverse cuts) was determined by 
consideration of the aspect ratio of the plate elements representing the 
tank skin. 

Th'i slosh baffles were constructed as a substructure with 208 
nodes; 104 were synchronous with the nodes on 8 of the LO^ tank shell 
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Figrure 2. Pictorial of model grid. 

transverse planes and 104 were unique nodes inside the tank on the 
same transverse planes. The slosh baffle structure was modeled using 
bar and plate elements. 

2. Intertank . The intertank is a setni-monocoque cylindrical 
structure with iflanges on each end and has a large box beam assembly 
(SRB beam) running transversely through the intertank (Fig. 3). 

The intertank and SRB beam were modeled as separate substruc- 
tures which were subsequently assembled. A remodeling of the SRB 
beam was necessary after test data demonstrated that the SRB beam 
modes were incorrectly predicted. The final model differed only in the 
representation of the cross section of the beam at the tank plane-of- 
symmetry which forms the half model boundary. 


4 



The elements used were plate elements for the skin and beam 
elements for flanges, longerons, and circumferential frames. Smearing 
techniques (grouping of stringers or ribs into a single finite element) 
were used to reduce the model size. 

The intertank nodal geometry was similar to the LO^ tank nodal 

geometry with grid points at the intersection of the 13 longitudinal planes 
with 6 transverse planes which represented the frames and flange rings. 

3. Support Ring . The support ring is a circular box beam 
designed to be rigid compared to the LC '2 tank /intertank assembly and 

to simulate the mass of an empty liquid hydrogen (LH^) tank. It mates 

with and attaches to the intertank flange ring provided for attachment 
to the LH^ tank. It was modeled as a substructure represented by 12 

prismatic beam elements. 

4. Airbag Suspension System . The suspension system consists 
of three clusters of airbags. One cluster spans the analysis plane of 
symmetry. The other two clusters are located at equal angular intervals 
from it. The suspension system was modeled as a substructure of eight 
nodes and six axial elements. Two of the nodes coincided with nodes 
common to the intertank and support ring. The axial elements were 
sized to provide stiffness characteristics of the airbag cluster. Those 
elements representing the airbag cluster which spanned the analysis 
plane-of-symmetry were halved to conform to the half-model consistently 
used. Since the number of airbags pressurized and contributing to the 
system support vtiried with the test configuration , a separate model was 
made for each configuration analyzed. 

5. Fluid Model . The fluid models were generated using a potential 
function approach with incompressibility being assumed. The method 
utilizes a surface grid. The tank wall nodes were defined to be synch- 
ronous with the structural model. In addition t these nodes a set of 
nodes was necessary to represent the liquid free surface. The surface 
was divided into elements bounded by the 13 longitudinal planes and 

3 concentric circles; the outermost circle is coincident with the tank wall. 
Two fluid models were necessary for each of the test configurations, one 
with symmetric boundary conditions applied and one with asymmetric 
boundary conditions applied. 


III. TEST CONFIGURATION ANALYSES 


Model solutions of the structural and fluid model of the test con- 
figuration were obtained from a MMC computer program, FORMA. The 
final solution degree-of-freedom set was obtained by a sequence of 
operations involving reduction, overlays, reduction and deletion. 
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This resulted in two half modal models of the LO^ tank /intertank /ring 

assembly supported on soft springs for each of the original four test 
conditions and the eight supplementary test conditions added during the 
testing. All of the modes up to 50 Hz were obtained. The resulting 
modal deflections and dynamic pressures were determined at the locations 
where test data sensors would be located. 

Hard copy plots of the resulting mode shapes were then mad s.jing 
an identical format as the test data plots. 


IV. TEST DATA 


A. Data Acquisition 

The test configurations investigated during the test program are 
described in Table 1. Wide band sine sweeps were made over the range 
of interest (3 to 50 Hz) and co-quad plots were made. Discrete fre- 
quencies and mode forms were identified from these plots. The MMC 
data evaluation team then chose those modes which they considered to be 
prime pogo oriented modes. The frequency of each significant mode was 
tuned until an acceptabk force- response /phase relation was achieved. 
Modal dwells and decay functions representative of the modes were then 
recorded on magnetic tape and on an on-line computer data storage unit. 


TABLE 1. LO 2 TANK TEST CONFIGURATIONS 


Configuration 

No. 

Liquid Level 
(in.) 

Cant Angle 
(deg) 

Ullage Pressure 
(psi) 

1 

487.0 

0 

3.3 

2 

320.7 

13 

3.3 

3 

162.0 

13 

3. 3 

4 

58.0 

13 

3.3 

5 

320.7 

0 

3.3 

6 

320.7 

0 

1.6 

7 

320.7 

0 

8. 3 

8 

320.7 

4 

3.3 

9 

384.5 

0 

3.3 

10 

384.5 

0 

8.0 

11 

218.0 

0 

3.3 

12 

162.0 

0 

3.3 


Modal dwell data were processed by the MSFC Computation Labora- 
tory into a tabular listing of the acceleration and phase response of each 
instrument and normalized mode shape plots of the test article for the 
longitudinal and lateral planes. 
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Modal decay data from selected instrumentation were processed by 
the MSFC test team into decay plots. 

Dynamic pressure responses in the LO^ tank aft dome were also 

recorded for well defined mode? . These data were normalii.efl to a local 
acceleration response and incorporated into the modal dwell tabuiar 
listing. 

Tables 2 through 5 present a summary of the test modal frequen- 
cies, their description and the percentage of critical damping for test 
configurations 1 through 4. 

Table 6 presents a summary of the dynamic pressure response of 
the LO,, tank aft dome for test configurations 1 through 4. 

Table 7 presents a summai y of the modal properties for test con- 
figurations 5 througl; 12. These are the supplementary tests that were 
performed to study the effect that fluid level, cant single, and/or ullage 
pressure combinations would impose on lightly damped modes. 

Table 8 presents a summary of the modal frequencies and damping 
variations due to cant angle/tluid level /ullage pressure changes on the 
first three symmetric bending and bulg'e modes for test configurations 
1 through 12. 


B. Data Evaluation 

A preliminary list of target modes, for which modal dwell and decay 
data were obtained , was compiled from pretest assessment of the analytical 
data for configurations 1 through 5. 

These data were reviewed by the MMC data evaluation team for 
validation of its modal character and acceptance. The primary criterion 
for its acceptance was the pha^e relation between the instruments and 
the primary force input. The tolerance established for phase deviation, 
the differential between the measured and theoretical phases, was 10° 
for all modes having high acceleration response. The secondary criterion 
was for high quality decay functions, a basic modal frequency absent of 
any harmonic modulation. The Fourier analysis method designed into the 
instrumentation computer system was employed as an aid to define any 
low response or potentially coupled modes. 

Tables 9 through 14 and Figures 4, 5, and 6 present summaries of 
the analytical and experimental data resulting from this program. Con- 
figuration 5 (Table 13) was introduced into the program as part of a 
supplementary test program required by information acquired in the test 
of configuration 1. An analysis of Ixie symmetric half-case was performed 
for inis configuration concurrent with the testing and provided additional 
test /analysis correlation data. 
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TABLE 2. TARGET MODE SUMMARY FOR 
TEST CONFIGURATION 1 




Test Modes | 

TSS No. 

Frequency 

(Hz) 

Da.nping 

(%) 

Description 


4.90 

1.6 

SYM SheU, M=i . N=2 


4.88 

1,6 

SYM Shell, M=l, N=2 


9.04 

1.2 

SYM First Pitch Bending 


5.72 

2.5 

SYM First Bulge 

028 

9.48 

.72 

ASM First Yaw Bending 

036 

12.76 

.3 

SYM Second Bulpe 

069 

18.95 

.^0 

SYM Third Bulge 


21.26 

.10 

SYM Fourth Pitch Bending 


19.29 

.18 

SYM Third Pitch Bending 

048 

16.63 

.28 

ASM Third Yaw Bending 

050 

21.55 

.12 

ASM Fifth Yaw Bending 

054 

23.83 

1.5 

SRB Beam, Z Bending 

056 

45.41 

1.0 

SRB Beam , X Bending 

058 

9.18 

.32 

SYM Shell, M=2, N=3 

059 

16.54 

.32 

SYM Lome Bending, Ogive Shell 

060 

14.08 1 

,17 

SYM Second Pitch Bending 

063 

25.75 

.43 

SYM Dome Bulge, Ogive Shell 

064 

19.68 

.51 

ASM Fourth Yaw Bending 


13.81 

.58 

ASM Second Yaw Bending 


9,75 

.i9 

SYM Shell, M=l, N=7 


TABLE 3. TARGET MODE SUMMARY FOR 
TEST CONFIG- 'lATION 2 




Test Modes | 

TSS No. 

Frequency 

(Hz) 

Damping 

(%) 

Description 

091 

6.64 

.60 

SYM Shell, M=l, N=2 


11.97 

.66 

SYM First Bulge 


13.02 

.44 

SYM First Pitch Bending 

098 

20.72 

.27 

SYM Third Pitch Bending 

099 

6.54 

.51 

ASM Shell. M=l, N=2 

105 

12.48 

.70 

ASM First Yaw Bending 

106 

19.82 

.76 

ASM Second Yaw Bending 

111 

19. 76 

.33 

SYM Second Pitch Bending 

114 

14.14 

1.2 

SYM Second Bulge 

125 

8.31 


SYM SheU, M=l, N=3 

126 

9.04 


SYM Shell, M=l, N=4 

127 

9.99 


SYM Shell, M=l, N=5 

128 

10.66 


SYM Shell, M=l, N=6 

129 , 

13.97 


SYM Dome Bulge, Ogive Shell 
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TABLE 4. TARGET MODE SUMMARY OF 
TEST CONFIGURATION 3 




Test Modes j 

TSS No. 

Frequency 

(Hz) 

Damping 

(%) 

Description 

149 

13.35 

.60 

SYM Dome Bulge 

150 

13.47 

.62 

SYM Second Bulge 

134 

20.04 

.50 

SYM Second Pitch Bending 

155 

14.91 

.90 

ASM First Yaw Bending 

156 

22.56 

.46 

SYM Distorted Dome Bulge, 

158 

16.54 

.60 

Ogive Bending 
ASM Second yaw Bending 

159 

19.36 

.90 

ASM Third Yaw Bending 

160 

25.37 

.57 

SYM Dome Shell, Ogive Bending 

168 

13.31 

.50 

SYM Shell, M=l, N=2 


TABLE 5. TARGET MODE SUMMARY OF 
TEST CONFIGURATION 4 




Test Modes [ 

TSS No. 

Frequency 

(Hz) 

Damping 

Ci) 

1 

1 

t 

Description | 

078 

15.67 

.23 

FYM ' 

j 079 

18.53 

.18 

Lienciin^* 

081 

;!5.5:! 

. 12 

SYM Shell and Pitch Bending j 


21.56 

.44 

ASM Second Yaw Bending I 

085 

18. '^8 

.41 

ASM Fi ■•si Y'aw .Bernting | 

087 

21. 12 

. 11 

SYM Second Pitch Bending ! 
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TABLE 6. DYNAMIC PRESSURE SUMMARY OF LO„ TANK AFT DOME FOR PRE-SELECTED 
















TAi5LE 7. TARGET MODE SUMMARY SUPPLEMENTARY TESTING 
rr£ST CONFIGURATION 5 THROUGH 12) 
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BLE ?. %iODAL FREQUENCE AND DAMPING VARIATIONS DUE TO CHANGE IN 
CANT ANGLE, FLUID-LEVEi-S AND ULLAGE PRESSURE (FIRST THREE TANK 
SYMMETRIC BENDING, AND DOME BULGE MODES) 
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TABLE 9. MODAL CORRELATION FOR TEST CONFIGURATION 
(487 in. DEPTH, 0° CANT-ANGLE) 
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TABLE 10. MODAL CORRELATION FOR TEST CONFIGURATION 
(320.7 in. DEPTH, 13° CANT-ANGLE) 
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11. MODAL CORRELATION FOR TEST CONFIGURATION 
(162 in. DEPTH, 13° CANT-ANGLE) 
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TABLE 13. MODAL CORRELATION FOR T -ST CONFIGURATION 
(320 in. DEPTH, 0® CANT-ANGLE) 
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Figure 4. Modes for configuration 





Figure 5, Modes for configuration 






Tables 9 through 13 present the test dwell data which were judged 
to have modal quafity. Table 14 presents the normalized modal pressure 
data for these modes. 

Part of this data is shown graphically in Figures 4 and 5. The 
plots also show shell modes identified from sweep data only. This pre- 
sentation is included to give insight into the trend of the test data 
relative to the analysis. The shell modes are identified by an M number 
and N number notation, where the M designates the number of half sine 
waves in the axia*. direction and the N designates the number of sine 
waves in the circumferential direction. The bending and bulge modes 
are identified as "first,” "second," etc., in order of ascending frequency. 

The data for the 54 test and corresponding analysis modes Isted 
in Tables 9 through 13 were subjected to careful detailed assessment in 
this test /analysis correlation. The assessment procedure started with 
examination 'f the analysis data for candidate target test modes. It 
continued tiirough the evaluation of the test data as describing valid 
modes. The establishment of corresponding test /analysis mode pairs 
was through detailed considoi-ation of the matching and the differing 
features of the mode shapes. The frequency match or mismatch was 
accepted after satisfaction of the shape criteria. 

Thirty-nine of the 54 mode p^rs were rated as having good or 
excellent correlation. The ratings were the result of subjective judg- 
ments of the correlation quality. Rigid criteria were not implied in 
arriving at the ratings as the relative importance of the various features 
concerned was different for different types of modes. Generally, the 
mode shape match of the lower dome was accorded the greatest weight. 
Thir was followed closely by frequency oori*espondence , then by mode 
shape match in the cylinder /ogive and finally by mode shape match in the 
int '•tank. 

Modal results from test configurations 1 through 4 identified 
several lightly damped modes within the 12 to 18 Hz frequency range. 

They were the symmetric bending and bulge modes. 

As a result of the SRB thrust oscillatory condition that was dis- 
covered at approximately 15 Hz during the SRM DM-1 and/or DM- 2 s'atic 
firings, supplementary test configurations were authorized for the purpose 
of studying the damping phenomena of these modes while varying the 
fluid level, cant angle, and ullage pressure parameters. 

Table 8 presents a summary of the modal frequency and damping 
properties of the first three symmetric bending and bulge modes for all 
test configurations (1 through 12). 

The trend established by those data indicates that modal frequency 
is a direct function of pressure and cant angle. It increases with higher 
pressure and with higher cant angle. Modal frequency, however, cannot 
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be directly related to fluid level variations due to the coupling effects 
of shell modes. No general trend can be established for modal decay . 

It can, however, be related to these variables on a per mode basis. 

To demonstrate this, part of this data is shown graphically in Figure 6. 

It is related to the second dome bulge mode and describes the modal 
decay as a function of these variables and its correlation within the 
thrust oscillatory time increment. The trends described by this mode 
during test were: 

a) Within the 0 to 320.7 in. level: 

1) It increased with higher fluid level for a constant pressure 
(3.3 psi) and cant angle (13^) test configuration. 

b) At the 320.7 in. level: 

1) It increased with higher ullage pressure for a constant 
fluid level (320.7 in.) and cant angle (0**) test configuration. 

2) It increased with higher cant angle for a fluid level (320.7 
in.) and pressure (3.3 psi) test configuration. 

c) Within the 320.7 to 487 in. level: 

1) It decreased with higher fluid level for a constant pressure 
(3. 3/8.0 psi) and cant angle (0°) test configuration. 

The trend previously described showed that the modal decay was 
low (0.13 percent) at liftoff as well as at the commencement of the poten- 
tial thrust oscillatory condition. It did, however, increase significantly 
during this interval, fium 0.13 percent to 0.78 percent at 3.3 psi or 0.14 
percent to 1.1 percent at 8.0 psi. 

Additional instrumentation was installed to the MVGVT test article 
to obtain an independent check of these damping values. The one data 
point that was obtained during the preparation of this report is shown 
in Figure 6 and correlates with the LO^ modal results. 

V. CONCLUSION 


It is concluded from the assessment of the LO^ modal teat and 

analysis data that the MMC hydroelastic finite element methodology is 
sufficiently accurate to duplicate the structural /fluid dynamic charac- 
teristics of a full scale flight external tank. The modeling technique 
used by MMC will accurately predict the primary "pogo" modes of the 
External Tank in all flight configurations. 
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